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FOREWORD 


This  report  was  prepared  by  the  Polymer  Branch  of  the  Nonmetallic  Materials 
Division.  This  work  was  initiated  under  Project  No.  7340,  “Nonmetallic  and  Composite 
Materials,”  Task  No.  734004,  “New  Organic  and  Inorganic  Polymers.”  It  was  adminis¬ 
tered  under  the  direction  of  the  AF  Materials  Laboratory,  Deputy  Cmdr/Research  & 
Engineering,  Aeronautical  Systems  Division,  with  Dr.  W.  E.  Gibbs  acting  as  project 
engineer. 

The  work  reported  here  concerns  the  solution  of  a  differential  equation  encountered 
when  treating  analytically  the  thermod)niamics  and  hydrodynamics  of  polymers  in  solution. 
The  present  work  represents  one  aspect  of  the  theoretical  and  experimental  program 
conducted  internally  on  this  topic.  This  program  has  as  its  general  objective  the  corre¬ 
lation  of  polymer  structure  with  physical  properties. 

This  report  covers  work  accomplished  from  July  1962  to  July  1963. 
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ABSTRACT 


The  following  Second  Order  and  First  Degree  Ordinary  Linear  Differential  Equation 
has  been  solved; 


X^y"  +X(  +  2t-4u)y'  +  (x^  +  2u-2t  +l)y 


(u  +  I)!  t !  2U  +  2 
(2t)!  ^ 


where  u  =  0,1.2,...®  and  t  =  0,1,2. ...u. 

The  general  solution  can  be  given  by  the  following  expression: 


y  = 


Hi 

-  I  (2 


!t! 


(2t)! 


+ 


dz  u 


[z  2 


+  u-t 


I 

•72 


(k 


I 


/zt-u-l 


54^  £jz 


>]}. 


where  z  =  x®,  and  k^,  kg  are  constants. 


This  technical  documentary  report  has  been  reviewed  and  is  approved. 


WILLIAM  E,  GIBBS 
Chief,  Polymer  Branch 
Nonmetallic  Materials  Division 
AF  Materials  Laboratory 
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DIFFERENTIAL  EQUATION  AND  ITS  SOLUTION 


The  following  set  of  differential  equations  (DE)  was  originated  in  a  treatment  of  hydro- 
dynamic  interaction  of  large  spheres  in  a  viscous  medium: 


L(X) 


X*y"  +  x(jX^+2t-4u)y'+(x^+2u  +  l-2t)y  = 


.2u+2 
(2t)!  ^ 


(I) 


where  u  =  and  t  =  0,1,2. ...u. 

1 

Introducing  a  variable  z,  where  x  =  z®. 


/ 

y 


2z  2 


dy 

dz 


and 


+  2 


dy 

dz 


,  equation  (1)  becomes: 


L(z)  =  4z 


2  d*  y 


+  z  z  +  4t  -8u  +2 


z  +  2u(2u  +  I  -2t  )  y  s 


(u+l)!t!  u  +  l 

e  —  V 


(2t)! 


(2) 


The  function  y  can  be  transformed  in  accordance  to  the  following  relationship: 


V  +  w  ,  where  w 


u!  t! 
(2t)!  ^ 


(3) 


Due  to  the  linearity  of  the  DE  under  consideration,  the  DE  of  y  can  be  expressed  as  a 
sum  of  differential  equations  of  v  and  w: 


/  ^  \  /  \  (u  +  l)!t!  „  +  | 

L(y(z))  =  l(v  (z))  +  l(w{z))  =  - ; —  z“ 

(2t)! 


(4) 


Since: 


2 

l(w(z))  =  4z^  +  zj^z+4t-8u+2j  +|^z+2u(2u-2t+l)Jw  , 


(5) 
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it  follows  that: 


l(w{z>)  =  [ 


4u(u-l)  +u(z  +  4t-8u  +  2)  +  z  +  2u(2u+l-2t)w 


,  .  .  .  .  (u  +  l)!t!  u  +  , 

(u+IJzW  =  - 7-r-  z 

(2t)! 


Therefore,  the  DE  of  v  is  homogeneous: 


l(v(2»)  =  42*  +  z[z  +4t-8u  +2] -^  +  j^z+2u{2u-2t+li]v  =  0  (7J 


DE  (7)  can  be  simplified  utilizing  the  following  transformation: 


V  =  z“  17 


where  tj  is  a  new  function. 


Since 


u  dTJ  .  u-l 

z  -rd-  +  u  z 


17  and 


14.  +  ia  + 

d  z*  dz 


/  11  ““2 

u(u-l)  z  77  , 


the  DE  of  the  function  rj  will  acquire  the  following  form: 


l(77(z>)  =  4z  — 4  +  (2+4t+2)  — -  +{u+l)^ 
dz  dz 


Introducing  a  variable  f  =  fz,  equation  (9)  becomes: 


=  C  ^  +  (C  +  t  +t)^  +{u+I)77  =  O 
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For  any  specified  value  of  u  and  t,  zero  or  positive,  a  constant 


k  =  (t  +  -^  )  -(u  +  l)  =  t-u-4- 


can  be  inserted  into  equation  (10): 


=  C  +  [C +k  +  (u +1)]  ^  +(u  +  l)-!7  =  0 


Differential  Equation  (11)  is  similar  to  the  equations  treated  by  A.  R.  Forsyth  and 
W.  Jacobsthal  (Lehrbuch  der  Differentialgleichungen.  Braunschweig  1912,  p.  759)  ref  (1). 
(See  also  reference  (2)).  An  analogous  treatment  of  equation  (11)  will  lead  to  a  special 
solution.  However,  an  additional  different  solution  will  also  be  found. 

Supposing  7]  ,  is  the  solution  of  DE  (11),  then  it  is  obvious: 

Uri 

'-u+|(’7u+|)  =  ^Vu+I  +  + '>]'9u  +  l  +  +  =  0 


The  differentiation  of  DE  (11 -a)  leads  to  a  new  equation: 


+  +  (u  +  2)i7'^ 


Hence: 


^  *-u+l^  ’I  u  +  l  ^ 


*"u+2^^u+2^’  '’7u+l  ”  "’7u  +  2  ond: 


‘-u+2^^u+2^  =  ^<+2  +  =  0  (12-0) 


Therefore,  knowing  any  solution  tj  of  DE  L  (rj  )  =  O  with  integer  S<u+1,  one  can  find 

s  s  s 

a  corresponding  function  by  consecutive  (u+l)-S  differentiations  of  and  such 

rj  .  is  a  solution  of  the  DE  L  ,  (tj  ,)  =  O. 

'u+1  u+1  '■  'u+1' 

Supposing  now  u  =  O,  the  corresponding  DE  will  be: 


L,(77,)  =  +  + ']  '^l'+  ’ll 
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The  differentiation  of  this  particular  case  will  lead  to  the  following  relationship: 

- =  C^,'"+[C+lc  +  2]i7"+2iy'  =  LgCijg*  = 

=  i-n"  +  =  0. 


where  %  =  tj/. 

It  is  sufficient,  therefore,  to  solve  DE  (13).  If  the  solution  %  is  known,  the  function 
can  be  found  by  u  differentiations  of  rji . 


It  is  important  to  mention  the  relationship  obtained  between  the  differential  equations 
of  type  will  be  valid  even  if  n  =  O.  Assuming  u  =  -1,  the  corresponding  DE  can  be 

written  as  follows: 


%  ■ 


To  prove  this  relationship  for  u  =  -1  it  is  sufficient  to  differentiate  the  DE 

-  + 'I’ll' =  L,  {77,1  =  O.  CI7I 


Equation  (16)  can  be  easily  solved.  Since: 


—  =  -(1+  |-)de, 
1,  « 


(18) 
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the  first  solution  of  tji  ,  ^  ,  will  be: 

=  C,  r'^  ,  (19) 


where  Ci  is  a  constant. 

To  prove  that  <f>  ^  also  solves  DE  (13)  and  <^”can  be  evaluated  and  inserted  into 
equation  (13). 

From  equation  (19)  one  can  obtain: 

J.I  .  /k  j//_  rk(k+l)  ,  2k  ..1, 

•Pi  -  J  +  '  )9,  and  =  [■  --f  +  '  j  ’/’i  ' 

Due  to  equation  (13),  if  ^  O  the  resulting  expression  will  be  valid: 

+  2k +C]-(|-  +l)(|  +  k +l)  + I  = 

=  +  2k  -It-I+I  :  0. 

i  i 


As  already  mentioned  (p,  is  not  the  only  solution  of  equation  (13).  The  second  solution 
<^2  can  be  obtained  if  Ci  of  equation  (19)  will  be  expressed  as  a  function  of  C  : 


-k 

=  C  e  c  (^ ) 


Since: 


^^2  _  ,  k  .  ,  X  j  .  £--k  d  C 

W  ■  ^  *  I?  ■ 


and: 


de^ 


_  rk(k+l)  2k  -1,  „/ k  .  ,x^-k  dC  .t-k-tliC 


-k -C  dc  ^  .-k  -e  ^ 


(22) 


{22-a) 


(22-b) 
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I  +  C']  = 


*“k 

Assuming  ^  O  one  can  obtain  from  equation  (23): 


Hence: 


^  -c 

de  ^  ^ 


ccC)  =  ca/e^Vde 


where  C,  is  a  constant. 


Finally: 


Since  no  other  solutions  exist,  the  general  expression  for  t?i,  is  a  linear  combination 
of  and  ^2  • 


%  =  C  ^  ^  ( C  +C  d^  ). 


Remembering 


^  ,  the  general  expression  for  ,  Is: 


^u+l  ^  dC' 
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As  k  =  t-u-|  this  solution  can  be  expressed  by  equation  (30): 


"7u+l 


dC 


e  '"{c^+c^  ^  ^  e^d 


(30) 


Utilizing  the  variable  z  =  4^  equation  (30)  becomes: 


^(ki+kg/z^  “"2  e^dz)] 


u+l  dz 


(31) 


where  and  kg  are  new  constants. 

Combining  equations  (31),  (8)  and  (3),  the  final  expression  for  y  appearing  in  DE  (2) 
is  given  by: 


y 


fu!  t! 
1(2 1)! 


e 


A. 


+  k 


-U--T  .-T 


dz 


(32) 


where  z  =  x^. 


TRANSFORMATION  OF  THE  FUNCTION  <f>z 


The  expression  for  given  by  equation  (30)  can  be  generally  written  as 


"7u+. 


(33) 


It  can  be  shown  that  the  power  of  C  appearing  in  the  integral  of  the  second  term  of 
equation  (33)  must  not  necessarily  include  u  and  t. 

Since  t  s  u,  then  O  ^  u-t=m. 

Substituting  m  into  the  expression  for  obtain: 


CgC 


•o  +  m 


'"e^dC 


(34) 
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The  integral  appearing  in  equation  (34)  can  be  transformed  as  follows: 


(in  +4) 


L_e«  +  + 


(m  +  j) 


df 


Repeating  this  procedure  m  times,  the  integral  under  the  consideration  becomes: 


u 


'■  s=o  (m  +  -r 


s=o  (m  +  -j)(m- 2)...(m  + 2 


+  4-S)  - 


(36) 


Utilizing  equation  (36)  the  function  can  be  expressed  by  the  following  formula: 


4>, 


,  m 

s=o 


■+m 


(m  +  ^)...(m  +  ^-S) 


<37) 


Since  s  ^  m  and  m  ^  u,  the  integer  s  appearing  in  equation  (37)  is  always  smaller  than 
u  with  the  exception  of  the  last  term  in  which  s  =  m  cannot  be  greater  than  u.  Therefore, 
differentiating  the  polynomial  of  equation  (37)  u  times,  all  terms  with  s  <  m  will  vanish 
and  the  term  with  s  =  m  can  be  zero  or  constant. 

Hence: 


fh 


A  +B 


8“ 


(38) 
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where  A  and  B  are  constants  and  A  =  O  if  m  <  u. 
If  m  =  u,  then 


A 


u: 


(39) 


Finally: 


(40) 
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